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Abstract 

This paper describes an algorithm for recovering a collection of linear cracks in a homo- 
geneous electrical conductor from boundary measurements of voltages induced by specified 
current fluxes. The technique is a variation of Newton’s method and is based on taking 
weighted averages of the boundary data. We also describe an apparatus that was con- 
structed specifically for generating laboratory data on which to test the algorithm. We 
apply the algorithm to a number of different test cases and discuss the results. 
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which is operated under National Aeronautics and Space Administration contracts NAS1-18605 and NAS1- 
19480; the third author is partially supported by National Science Foundation Grant DMS-9202042 and 
AFOSR Contract 89-NM-605. 


l 




1 Introduction 


The purpose of this presentation is to demonstrate the viability of using Electrical Impedance 
Tomography (EIT) for the reconstruction of multiple macroscopic cracks. Through the use 
of EIT one seeks to determine the interior conductivity properties of a specimen from mea- 
surements of electrostatic potentials and currents on the boundary. There are obvious appli- 
cations of this technology to many problems, for instance in medical tomography [5], [8] or 
nondestructive testing of mechanical components [10]. There are two very separate classes of 
practical reconstruction techniques: 1) general imaging techniques, and 2) techniques that 
use prior information or assumptions. The first approach produces a somewhat blurry image 
of the conductivity distribution inside the specimen (somewhat akin to an ultrasound-scan) 
from which an “educated” user may then draw conclusions of a more specific nature (c/. 
[5], [8], [10], [18]). The second approach uses a restricted model for the interior conductivity 
distribution and seeks to determine more specific features within this model which are con- 
sistent with the measured data (c/. [13], [15], [16]). This may be done probabilistically, in 
a way reminiscent of methods for image reconstruction known as “maximum entropy meth- 
ods” or deterministically by seeking to fit a relatively small number of parameters. In this 
presentation we shall concentrate on a very specific method of this last category. It seeks 
to determine a finite number of linear cracks consistent with the boundary measurements. 
In particular we shall test this algorithm on data which has been collected from laboratory 
experiments. 

The mathematical results which insure that a finite number of cracks can be reconstructed 
from a finite set of boundary voltages and boundary currents are found in [6] and [11]. In 
these papers it is proven that n C 2 cracks inside a two dimensional domain may be determined 
from n + 1 pairs of boundary voltages and boundary currents. These results have recently 
been extended to show that measurement of two pairs of boundary voltages and boundary 
currents suffice to determine any number of cracks [2], [9]. It is very easy to see that one 
measurement will not suffice to determine even a single crack [11]. 

The computational algorithm we shall employ has largely been developed in [7]. It is the 
natural extension of an algorithm which was developed in [17] for the case of a single crack. 
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We have modified the algorithm from [7] in one important respect, by making the selection of 
the initial guess automatic. For this purpose we have implemented some formulae which were 
suggested in [4] and which directly reconstruct the line on which a single linear crack must 
lie in order to be consistent with a single pair of boundary voltages and boundary currents. 
These formulae may occasionally be degenerate, but the use of additional measurement pairs 
can eliminate this degeneracy. 

An outline of this presentation is as follows. In section 2 we introduce the mathematical 
model used for simulating cracks and we give a brief discussion of our computational recon- 
struction method. In section 3 we present the formulae derived by Andrieux and Ben Abda 
and we briefly describe the selection of the initial crack. Section 4 contains a description of 
the experimental equipment used to collect the actual data for our reconstructions, and in 
particular we describe the few modifications that were made when compared to the equip- 
ment used earlier to collect data for the reconstruction of a single crack [17]. In the fifth and 
final section we provide a number of examples of reconstructions of multiple cracks based on 
the collected experimental data. 


2 The Mathematical Model and the Reconstruction 
Algorithm 

A single crack inside a two dimensional conductor is commonly modeled as a perfectly 

insulating curve < 7 . With a background conductivity 0 < 70 < 7 (z) < 71 and a finite 

collection of cracks S = UJLjtfjfc, the steady state conductance equations thus read 

V * ( 7 Vu) = 0 in fi\E, (2.1) 

dv 

7 ^ = 0 on E, 

with appropriate boundary conditions on dfl, e.g., 

v = <j> on dfl. (2-2) 

The field v is normal to E. The function v represents the potential induced in fi. We assume 
that C 1R 2 is simply connected, i.e., has no holes, and so the entire boundary dQ is 
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accessible from the “outside”. The domain corresponding to our experimental data is a 
disk. Let u denote the “7-harmonic” conjugate to v. It is related to v by the formula 

(Vu) 1 = jVv in fl \ E, (2-3) 

where ± indicates counterclockwise rotation by 7r/2. Note that the set fl \ E is not simply 
connected; the existence of a “7-harmonic” conjugate u owes not only to the fact that 
V • (jVv) = 0 in fl \ E, but also to the fact that 7^ = 0 on E. For a particular set of 
constants q t, k = 1, ... , n, the function u solves the problem 

V ' (7 -1 Vu) = 0 in n\E, (2.4) 

u = cjt on 07, k = 1 , ... ,n 

with 

= on an. (2.5) 

ov Os 

Here s denotes the counterclockwise tangent direction on dfl and v denotes the outward 
normal on dCl. For these particular constants, finding a solution to (2.4)-(2.5) is thus equiv- 
alent to finding a solution to (2.1)-(2.2) . The constants c k may (up to a common additive 
constant) be characterized in several equivalent ways (see [7]). The characterization we shall 
use here takes the form of the system of n equations 

[ [7 -1 ^] ds = 0, 1 < k<n. (2.6) 

J <jk uu 

Here [7 -1 §^] = 7 _1 ^r — 7 -1 denotes the jump in the normal flux across the curve a k . 2 
The system (2.6) should be viewed as a supplement to the boundary value problem (2.4)- 
(2.5). The solution u has a physical interpretation in its own right - it is the potential 
generated by the boundary current xp for a medium with background conductivity 7 _1 and a 
set of perfectly conducting cracks 07, k = 1, ... ,n. As follows immediately from the relation 
(2.3), one passes between the boundary data corresponding to v and that corresponding to 
u by differentiating the Dirichlet-data along the boundary, integrating the Neumann-data 
along the boundary and interchanging the roles of the two. Except when we explicitly say so 
2 The expression denotes the limit of the derivative (in the direction v) as one approaches rr k from 

the side to which v points. 4^- denotes the limit as one approaches cr k from the opposite side. 
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we shall always work with cracks that are perfectly conducting in the sense described above. 
Let P\ , ... , Pm and Q\, ... , Qm be 2 M points on dVl. For the crack reconstruction we utilize 
solutions corresponding to the two-electrode currents tf>j = 6 p } — j = 1 , . . . , M, 

0 in fl \ E, (2-7) 

<£> on <7*, k = 1 , ... , n, 

7 _1 = ^P, ~ &Q, on dCl, (2.8) 

the constants cjj.^ being determined through the additional equations (2.6 ). The reconstruc- 
tion problem may now be stated explicitly as follows: 

We seek to reconstruct the collection of cracks E — U£ =1 <77 from knowledge of the boundary 
voltage data {wj ianl^i corresponding to the prescribed two-electrode currents 7 -1 = 6 p i — 

j — 1 > ••• > M . 

It was shown in [6] that if we take M = n + 1, Qi = Q 2 ... = <? n +x = Po and we take 
Pj, 1 < j < n + 1 to be mutually different and different from Po, then the boundary voltage 
measurements corresponding to the resulting n -f 1 fixed two-electrode currents 8 P, ~ tPoi 
1 ^ j ^ n + 1 suffice to uniquely identify a collection of n (or fewer) cracks. By a clever 
extension of this analysis it has recently been shown that two fixed two-electrode currents 
suffice to uniquely identify any number of cracks, [2], [9]. It is very easy to see (cf. [11]) that 
measurements corresponding to a single boundary current do not suffice to determine even 
a single crack. In the papers [6], [11] it was required that the reference conductivity 7 be real 
analytic. This requirement has recently been considerably relaxed; the identifiability result 
holds even if 7 is only Holder continuous [9]. 

A certain amount of knowledge is available concerning continuous dependence. In the 
case in which the background conductivity is constant and there is only one crack present 
it has been shown (cf. [1]) that if the boundary voltage data (on some open subset of df!) 
deviate by e, then the crack locations differ by at most C[log(| log e|)] -1 ^ 4 . If the single 
crack is linear it has been shown that the crack location depends Lipschitz continuously on 
the boundary data, [3]. This last result significantly extends a calculation found in [11], 
showing that transverse translations and rotations of a single linear crack depend Lipschitz 


V • (l-'Vf,) 


u, = 


with 


. flu 
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continuously on the measured data. 

As mentioned before, we shall in the present reconstructions always try to fit the data 
by means of linear cracks. We shall also assume that the background conductivity may be 
modeled by a constant, e.g 7 = 1. As was the case in [17], we base the reconstruction of 
the cracks on the values of a relatively small number of functionals (as opposed to all the 
boundary measurements). In [17] we used 4 functionals for the reconstruction of a single 
linear crack, here we use 4 n functionals for the reconstruction of n cracks. We now give a 
brief description of these functionals. For more details we refer the reader to [7], [17]. 

Let F denote the vector-valued function 

F(E, xp, w) = (F(E, rp, u; (1) ), F( E, ip, ™ (2) ), F(E, xp, u> (3) ), F( E, xp, u> (4) ))‘, 


where F(Y,,xp,w) is given by 

F(E, xp, w) = jf ^ u(E, ds. (2.9) 

The functions 1 < i < 4, will be taken as particular solutions of 

A u» = 0 in ZR 2 \ E. 


The function u = u(E, xp) is the solution to (2.4)-(2.5) with 7 = 1. In order to make u unique 
we may for instance impose the requirement 

[ uds = 0 . (2.10) 

Jd n 

The appropriate selection of boundary currents xp and test functions w = (u^\ u/ 2 \ xv^.w^) 1 
is very important and was discussed in detail in [7]. For our reconstruction algorithm we 
always choose xp in the form of a two-electrode current. We take one xp and one w corre- 
sponding to each crack whenever we want to emphasize this correspondence we use the 
notation 

xpa k and = (u^ , w ¥ k » w ? k > w ^ k )*• 


We shall always take w so that 


L 


da dv 


ds = 0, and 


f 

dxo^ 

L 

dv 


ds = 0 Vcr G S. 


( 2 . 11 ) 
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Because of the first identity in (2.11) the function F is unchanged by the addition of a 
constant to tt, and we can therefore work with any other normalization in place of (2.10). 
The data for our reconstruction consist of measured boundary potentials. In practice we can 
of course only measure the values of these potentials at a finite (fairly moderate) number of 
points. We denote by g(xp) the voltage data corresponding to the boundary current ip, and 
we define a corresponding vector-valued function 

f (0,w) = (f(^w {1) ),f(ip,w {2) ),f(ip,w i3) ),f(ip,w i4) )y, 

where f(ip, w) is given by 

fW> w ) = j da 9 W^;ds, ( 2 . 12 ) 

and where w W are the same functions as before. Our algorithm seeks a solution S = {crfc}£ =1 
to the 4 n equations 

F(S,^ fc ,w <r J = f (^,w„J, l < k < n. (2.13) 

We do not use information about the full set of measured boundary voltages for the recon- 
struction; we only use information about the values of these particular functionals. This is 
in contrast to the reconstructions from experimental data which we presented in [13], and 
where we used all the measured boundary voltages and a least squares approach. Based 
on extensive experimentation it is our experience that almost all the relevant information is 
contained in these functionals. Admission of the extra data, e.g. after convergence of the 
algorithm based on the functionals, does not generally seem to improve the reconstructions. 
In several examples in which we tried the least squares approach from the start it actually 
impeded the convergence process, which we think is due to the presence of local minima. 

The construction of the test functions u>h) i s quite simple to describe. If we consider a 
single linear crack, and use a coordinate system so that the crack lies on the x-axis (the line 
of reals) then the functions u/ 1 ) and are given by 

= Im[z\, w™ = Im[z 2 ] . 

The functions w ^ and have square-root singularities at the endpoints of the crack; if 
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the crack lies between the origin and the point (A,0) they have the form 

M* - A)\A(z - A)], Re(z) > A u 

-fie[(z - A )yjz(z - A)], He(z) < A , 

w <*) = I Re{ ^ Z ~ A)1, ^ > ’ (2.15) 

-Re[^Jz{z — A)], J2e(z) < | . 

Intuitively w ^ and detect transverse translations and rotations of the crack, whereas 
u;( 3 ) and te (4) detect translations along the direction of the crack and variations in length. 
It is interesting to note that in the limit as A approaches zero the information extracted by 
these four functions approach the first four nontrivial Fourier modes of the boundary data. 
We refer the reader to [7] for more details. 

The selection of the applied currents is also fairly simple to describe: for a single crack 
we select the location of the two electrodes in a way which maximizes the sensitivity of the 
four functionals with respect to rotations and transverse translations of the crack. We use 
an iterative approach to solve (2.13), in this case for the single crack a. At each step new 
electrode locations are found by maximizing two of the diagonal elements of the Jacobian 
of D„ F with respect to the boundary electrode locations (taking the fixed test functions 
described above). Measured data corresponding to these electrode locations are then used 
to find an updated crack. The selection procedure is described in detail in [7], [17]. When the 
domain is a disk, and the crack is of moderate size and sits away from the boundary, then 
the most sensitive electrode locations defined above turn out (very nearly) to be the points 
that arise as the intersection of the line on which the crack lies and the domain boundary 
(remember, the crack is perfectly conducting). Our selection of electrodes does require the 
solution of an auxiliary boundary value problem. However, as explained later, this is done 
very efficiently through an integral equation formulation. For the reconstruction of more 
than one crack we select a boundary current corresponding to each crack, as if the other 
cracks were not present. This process is described in detail in [7]. 

The 4n x 4n system (2.13) is solved by means of a Newton method. As explained above, 
the test functions and the applied currents change with each step of the iteration. The 
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update for each step of the Newton method is constrained, along the lines of the classical 
Levenberg-Marquardt approach [14]. 

The boundary value problems which must necessarily be solved in order to evaluate F 
(and find the “optimally sensitive” electrode locations) are formulated in terms of boundary 
integral equations. The system of boundary integral equations is a mixture of first and 
second kind equations. There is a first kind contribution corresponding to the functions 
that describe the normal derivative “jumps” on each crack. However, by working with the 
boundary condition u =constant on each crack, we avoid the problem associated with non- 
integrable kernels. We di&cretize the integral equations by Nystrom’s method, as completely 
described in [7]. The kernels of the first kind contributions are singular enough that we 
did not find it necessary to' regularize the equations. Also, we note that we do not really 
numerically compute the it|ao necessary to evaluate F, rather we compute (u — u 0 ) |aa, where 
uq is the potential generated by the exact same electrode locations and currents as u, but 
without any cracks present. The computation of the smooth function (u — tto)|dn is a much 
better posed problem than the computation of u|an- Because the domain 0, is a disk, the 
function uo has a very simple explicit expression. Whenever necessary we add its values to 
the computed differences. 


3 Apparatus and Data Collection 

In this section we give a brief description of the apparatus with which we collected the 
experimental data used to test the crack reconstruction algorithm. Most of these details can 
also be found in [13]; nonetheless, for the convenience of the reader we will recount them 
here. 

The reconstruction algorithm is designed for a two dimensional region with a constant 
background conductivity. An experimental equivalent is given by the tank illustrated in 
Figure 1. The tank is cylindrical and constructed of 0.25 inch thick plexiglass, with an inside 
diameter of 10.5 inches and a height of 23.5 inches. The tank is filled with distilled water 
to which a small amount of ordinary tap water is added, which gives the solution a modest 
electrical conductivity. Twelve evenly spaced copper electrodes run vertically along the full 
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Figure 1: Tank used for data collection. 

height of the tank. The copper electrodes are 0.25 inches wide. Since the entire apparatus is 
uniform in the vertical direction, and since the electrodes are essentially perfectly conducting 
(in relation to the interior solution), the tank provides a reasonable approximation to a two- 
dimensional electrical conduction problem. To introduce perfectly conducting “cracks we 
insert strips of sheet metal with widths ranging from 2 to 6 inches vertically into the tank, 
again preserving the uniformity of the conduction problem in the vertical direction. A grid 
is placed under the tank (which has a plexiglass bottom) so that the location of the cracks 
can be easily recorded. 

The overall structure of the data collection apparatus is as follows. To each electrode 
is attached a binding post on the top rim of the tank. These binding posts are connected, 
through a bank of relay cards, to a voltage source and a pair of digital multimeters. One 
multimeter serves to measure the electrical current supplied to the tank when a specified pair 
of electrodes is active. The other multimeter measures the induced voltage drop between any 
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given pair of electrodes. By altering the electrical connections via the relay cards, any pair 
of electrodes can act as the source and sink for the electrical current and the voltmeter can 
measure the resulting potential across any two electrodes. The entire apparatus is connected 
to a table-top computer (PC) which controls the relay cards, allowing software switching of 
the active electrodes and voltage measurements. The PC and software also communicate 
with the multimeters, collect the resulting current and voltage information, and present it 
in a convenient format. A simple schematic of the entire apparatus is given in Figure 2. 



Figure 2: Schematic of data collection system. 


As mentioned in [13], while the problem is modeled as a steady-state or DC conduction 
problem, in practice one uses an AC current. This is to avoid any electrochemical plating 
of the electrodes with the impurities present in the water, which might alter the behavior of 
the electrodes, particularly the contact resistance between the electrodes and water. We use 
a sinusoidal source with a frequency of 10 kHz, which should be low enough to be considered 
steady-state so that no phase changing capacitance or inductance effects will be significant. 
In actuality, a few experiments were conducted with frequencies from 10 Hz to 1 kHz without 
any cracks in the tank, and in this range no frequency dependent effects were observed. 

The data collection software is written in BASIC. As discussed in the previous section, 
at the heart of the crack recovery algorithm is the idea of iteratively adapting the active 
electrode locations for greater sensitivity as the algorithm progresses and more precise infor- 
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mation about the actual crack location (s) becomes available. However, since the reconstruc- 
tion algorithm could not immediately be ported to the PC which collects the data, it was 
necessary to collect all possible data a priori. When the reconstruction algorithm was later 
run using this data set (on a separate computer), it simply extracted the required voltage 
data corresponding to the desired active electrode locations. One can easily check that for 
a tank with 12 electrodes there are only 11 different linearly independent combinations of 
input and output electrodes; the solution corresponding to any other current pattern can 
be represented as a linear combination of the solutions corresponding to these basic cur- 
rent patterns. In its dipole current pattern mode the software collects data for the current 
patterns = ^i+i ~ i = 1, . . . , 11, that is, a current is introduced at electrode j + 1 
and withdrawn at electrode j. Whenever the recovery algorithm needs a different current 
pattern during its execution it simply takes an appropriate linear combination of the data 
corresponding to these 1 1 basic patterns. 

For the data sets which will be presented, typical values for the input current were 4 to 
7 mA (rms). Typical voltage drops between electrodes range from zero to 7 volts rms, with 
the largest drops between the active input and output electrode. The previous version of 
the software, used to collect the data in [13], recorded the input current to only 2 significant 
figures. In the course of taking measurements for this paper it was observed that even the 
third significant figure for the current was quite consistent and stable from measurement 
to measurement. The software was thus modified to record three significant figures for 
the input current. As before, voltages are recorded to four figures. An application of the 
crack reconstruction algorithm to both types of data (two versus three significant figures for 
current measurement) indicate that the latter is considerably more accurate. It appears that 
for this experimental setup, the crack manifests much information concerning its presence 
and location in this third significant digit, particularly for multiple crack problems. 

Before proceeding with the reconstruction algorithm one must also know k, the back- 
ground conductivity of the tank. To obtain this information we applied the current patterns 
Vb, • • • , 0 ii described above to the tank with no cracks and measured the corresponding in- 
duced voltages. We then used the same current patterns in our computational model and 
computed the corresponding theoretical voltage. These voltages scale linearly with and so 
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one can easily compute that value for k which gives the best fit of the computational model 
to the measured voltages, in the sense of mean square error. Doing this yielded an average 
resistivity of 456.6 ohms, or a conductivity of 0.00219 mhos. This value varied by only 1 
percent over the course of the tests (approximately 30 different crack configurations over 2 
days). This is the value of k which was used in all crack reconstructions. In practice we 
scaled the measured data by k so that it formally corresponds to a background conductivity 
of one. All our subsequent computations were then done with a conductivity of 1 . 

One final remark is in order. The electrodes on the actual tank have a necessarily non zero 
width, while the simple theoretical model we use assumes point electrodes. This assumption 
leads to logarithmic singularities in the theoretical voltage at the current input and output 
electrodes. Moreover, in the* actual experimental data any contact resistance between the 
electrode surface and the water leads to inaccuracies in the voltage measurements at the 
active electrodes, since Ohm’s law implies that if a current is flowing through the electrode 
and a resistance is present, there must be a voltage drop across the electrode. Thus the 
voltage on the tank side of the active electrode (the voltage required by the algorithm) may 
differ from the voltage measured by the multimeter which is on the other “side” of the contact 
resistance. This is not a problem for voltage measurement s 1 at non-active electrodes; these 
electrodes have essentially no current flowing (the multimeter presents an input impedance 
of > 10 9 ohms), hence no voltage drop. We thus omit from consideration the voltages at 
active electrodes for both the experimental and computational data. One merely adjusts 
the numerical integration rule for computing the functionals F(£,^,u;^) and /(^,u>W) to 
account for the gap. 

While a more accurate model incorporating a non-zero width electrode could certainly 
be derived and the algorithm adapted to deal with such a model, the agreement between 
the measured data and theoretical computations was so good that there seemed little to be 
gained by refining the model, at least for cracks which are not extremely close to an active 
electrode. 
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4 A Simple Formula for the Line of a Single Crack 

In [4] two very simple formulae are given for reconstructing the line on which a single linear, 
perfectly insulating crack lies. These formulae are in terms of a single voltage/current pair 
measured on the boundary of 0. The formulae are restricted to the reconstruction of a 
single line. As has already been observed, a single crack (even a linear one) is not uniquely 
identifiable from a single measurement. As far as these explicit formulae are concerned this 
expresses itself in terms of the possibility of the vanishing of a certain denominator (more 
about this later). The formulae do not require the solution of any boundary value problems, 
but only the evaluation of certain weighted integrals of the voltage and current data on the 
boundary of the region. These formulae merely locate the line on which the crack lies, they 
do not tell us where on the line the crack is. It is possible to determine the exact location of 
the crack by evaluating an infinite number of similar integrals and reconstructing the support 
of a function from its Fourier series (see [4]). We have not implemented this latter step as 
part of our algorithm, but we have incorporated the simple formulae for the reconstruction 
of a single line as part of the initialization step of our crack reconstruction algorithm. Once 
the line is found, the initial guess is chosen to be centered on that segment of the line which 
lies inside the domain (a disk). In the reconstructions we present later the initial guess is, 
rather arbitrarily, taken to be of length 0.1 x the radius of the disk. 

We now briefly review the formulae from [4]. Suppose that cr is a perfectly insulating 
linear crack in and u(x) is a function on fl \ S C IR 2 which satisfies 

Av = 0 in ft \ a 

dv 

— = 0 on cr , 

dv 

where v is a unit normal vector field on cr. Note that this field is constant ( v = Vo ) since 
a is linear. Based on the boundary values of v and on d£l we can then recover the line 
vo • x — c on which the crack cr lies as follows. Let <j>i{x) denote any of the two coordinate 
functions <f>i(x) = X{. Note that <f>i are both harmonic. Define the quantities i = 1,2, by 

* =/«,(” * ( * )= *- (41) 

These quantities depend only on boundary values of v and its normal derivative, and can 
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be easily computed from the known data, e.g., when is specified (a current applied) and 
v is measured on dCl, or vice-versa. Let [v] denote the jump in u across <7; [v] is defined as 
— u_, where v + denotes the values of v as one approaches a from the side to which the 
normal vector Vo points, and u_ denotes the values as one approaches from the opposite side. 
A simple application of the divergence theorem shows that 


, ' = l [v] ib ds ' 

since vanishes on <r. Now note that ^ = v 0i i on <7, since fa = x^. Thus 


/, 



(4.2) 


If we assume ^ 0 it then follows that 

iL^ ds ) = 


and 

i/o,,' = ±Ii/yfl } + /| • 

This gives us a normal vector to the line on which a lies. To be precise let us choose vq 
corresponding to the + sign, i.e., 

vo,i = U /yflf+Tl ■ (4.3) 


From (4.2) it follows that 


for this choice of i/q, and therefore 



ds > 0 


J [v] ds = \J 7j 2 + I\ . 

We now recover the constant c, of the equation i/ 0 • x = c for the line on which a lies. Let us 
rotate our coordinate system so that the normal vector vo to a has coordinates (0, 1). Since 
we know j/ 0 this rotation is known; the line on which a lies is still u 0 • x = c, or simply x 2 = c. 
Let <j>(x) = \{x\ — xj) in this new coordinate system. The function 4>{x) is again harmonic. 
Define by 


86 , dv 



Application of the divergence theorem shows that 



Since V<f> = (— x\,xi) and since vq = (0, 1), this becomes 



or 



so that 

c = /*/(jf [v] <!») = /*. /•//? + /?, 

again assuming that the denominator does not vanish. This formula in combination with 
the formulae (4.1) and (4.3) gives a simple way of determining the line on which a single 
linear, perfectly insulating crack lies. 

We cannot generally take for granted that f^v] ds is non-zero, given a single fixed bound- 
ary current on 5ft. This is certainly the generic case, but it is not difficult to see that a 
single applied boundary current can not guarantee that this is true for all a . This is consis- 
tent with the general result that one needs two applied currents and corresponding voltage 
measurements even to find a single crack. If the applied currents are all two-electrode cur- 
rents it is very easy to see that a finite number will suffice in order to obtain I\ + /| ^ 0 
for any interior crack for at least one of these currents. By a more detailed analysis it is 
apparently possible to show that any two different two-electrode currents always suffice [2]. 

We note that the choices for the three auxiliary functions used for reconstructing the line 
on which the crack lies are harmonic polynomials of degree < 2, quite similar to the weight 
functions used in our crack reconstruction algorithm. In the limit (as crack length goes to 
zero) our four weight functions u/ 1 ) through span the space of harmonic polynomials of 
degree < 2, modulo constants. This is also the reason why in the case of a disk the limiting 
information contained in the functionals is equivalent to the first 4 nontrivial Fourier modes, 
as noted earlier. 

Finally let us note that it is very easy to write the quantities /,■ and I $ in terms of 
data corresponding to the harmonic conjugate to v, i.e., in terms of data corresponding to a 
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perfectly conducting a. Given a smooth function / on dil which satisfies f gn fds = 0, let f s f 
denote a counterclockwise integral of / along dCl, i.e., a function whose counterclockwise 
tangential derivative equals /. From the relation (2.3) with 7 = 1 we now immediately 
obtain 


where ^ denotes the counterclockwise tangential derivative on dtl. Similar formulae for /,- 
are obtained by replacing (f> with fa. We use these formulae directly instead of making any 
transformations to our measured boundary data. 


5 Reconstruction Experiments 

We have performed numerous experiments with the equipment described in section 3 followed 
by reconstructions using the algorithm described in section 2. In this section we discuss the 
results in three different scenarios, that we find to be representative. In the first scenario 
two sheet metal strips were inserted in the tank to simulate perfectly conducting cracks. We 
varied the locations of the two metal strips inside the tank in order to give some idea of 
how the information content of the measured data deteriorate as the cracks lie deeper and 
deeper inside the domain. In the second scenario three sheet metal strips were inserted to 
demonstrate that the data do indeed permit one to distinguish more than one or two cracks. 
We have not at the moment experimented with a number of cracks larger than three, partially 
because of the practical difficulties associated with the accurate insertion of a large number 
of cracks, but also because we feel it is likely that the currently available accuracy of the data 
(and the mathematical model) will make such reconstructions extremely crude. In the third 
scenario we used a single metal strip, but this time bent so that its cross section forms an 
L-shape. This experiment demonstrates that the collected data have sufficient information 
content to distinguish a shape which differs from a line segment, and it also demonstrates the 
flexibility of the reconstruction algorithm to approximate such a shape. In our opinion the 
results presented here clearly demonstrate the viability of impedance imaging as a method 
to reconstruct a moderate number of macroscopic cracks from experimental data. 
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Let us give a brief description of the figures. Each figure depicts two disks to illustrate 
one iteration of the reconstruction algorithm. The location of the reconstructed crack(s) are 
indicated by solid lines. Those in the disk to the left are the locations prior to the particular 
iteration, and those in the disk on the right are the locations after the iteration is complete. 
The dashed lines in the disk to the right indicate the locations of the “true” cracks as read 
off from the mesh at the bottom of the tank. The “optimally sensitive” electrode locations 
that were used for this iteration are indicated on the boundary of the disk to the left. We 
select one pair of “optimally sensitive” electrodes corresponding to each crack, according to 
the strategy briefly described earlier in section 2; n cracks will thus in general correspond to 
2n active electrode locations. We do permit some of the active electrode locations to coincide 
- though we always insist that the n prescribed two- electrode boundary currents be linearly 
independent (so there must be at least n + 1 distinct active electrodes). This strategy is 
different from that of [7], where we also selected one pair of electrodes corresponding to 
each crack but fixed one electrode location to be shared by all the electrode pairs; in those 
reconstructions n cracks corresponded to exactly n + 1 active electrode locations. We found 
that in the presence of the “noisy” experimental data the reconstructions were improved 
by not forcing one common electrode. Although we have not indicated on the figures which 
electrode pair corresponds to which crack, it is almost evident, given the observation we made 
earlier concerning a single perfectly conducting crack. The optimally sensitive electrodes are 
located near the intersection of the boundary and the line on which the crack lies. Figure 
3a shows the first iteration using experimental data corresponding to two cracks, but in this 
case the algorithm seeks to reconstruct only a single crack. The initial guess shown in the 
disk to the left lies on the line generated by the formulae of Andrieux and Ben Abda. We 
have chosen the initial guess to be of length 0.1 xthe radius of the disk, and to be centered on 
that segment of the line which lies inside the disk. Figure 3b shows the fifth iteration of the 
algorithm when starting with the initial guess shown in Figure 3a. At this point the algorithm 
has converged in the sense that the maximum of all four components of G is less than 10 10 
in absolute value. We now take this crack and divide it into two by omitting a central piece, 
1/10 of its size. The resulting two cracks are given as initial guess to our algorithm for the 
reconstruction of two cracks. After 83 iterations the reconstructed cracks have converged 
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The first iteration reconstructing one crack. The “true” crack locations are shown in dashed 
lines in the disk to the right. The currently active electrode locations are indicated on the 
circle to the left. 

to the two cracks shown in the right disk of Figure 3c. We have run several examples with 
data originating from two cracks. As one would expect the reconstructed cracks deviate 
more from the “true” cracks when the “true” cracks are closer to the center of the disk. The 
above reconstruction is quite representative of cracks located about 1/3 of the way towards 
the center. Figures 4a and 4b show the iterations at which convergence is reached when 
reconstructing one and two cracks respectively, given data that correspond to two “true” 
cracks located nearly halfway towards the center. Figures 5a and 5b show reconstructions 
obtained from data corresponding to two “true” cracks located approximately 3/4 of the way 
towards the center. Figure 5a shows the fifth iteration reconstructing only a single crack, 
at which point this process has converged. We divide the resulting crack into two pieces by 
cutting out a central piece 1/10 of its size, and feed these two cracks as initial guess to our 
algorithm reconstructing two cracks. Figure 5b shows the 90th iteration of the two-crack 
reconstruction process. There is no convergence, nor does any converged state appear to 
be reached at a later iteration. The collected data is apparently not accurate enough to 
distinguish these two cracks from a single crack near the center. We now proceed with a 
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The 5th iteration reconstructing one crack. The collected data are the same as in Figure 3a. 
The “true” crack locations are shown in dashed lines in the disk to the right. The currently 
active electrode locations are indicated on the circle to the left. 

scenario involving three “true” cracks. We start by reconstructing one crack which, in the 
sense of our functionals, is consistent with the data. The initial guess is generated by means 
of the Andrieux-Ben Abda formulae, as described earlier. We take the converged crack as 
shown in Figure 6a (after 6 iterations), divide it into two and feed the resulting cracks as 
an initial guess to our algorithm for the reconstruction of two cracks. After convergence (in 
this case after another 33 iterations, as shown in Figure 6b) we divide the largest of the 
reconstructed cracks into two pieces, and feed the 3 resulting cracks as an initial guess to our 
algorithm for the reconstruction of three cracks. Figure 6c shows the results of the three- 
crack reconstruction after 90 iterations. Convergence has not been attained, but notice that 
the reconstructed cracks none the less give a reasonable prediction of the “true” locations. 
Also notice that there are only 5 active electrodes on the left circle in Figure 6c. This is 
therefore an example where two active electrode locations coincide, as discussed earlier. 

Finally, in Figures 7a and 7b we show the converged reconstructions using one and two 
cracks respectively and data which come from insertion of a single strip of metal which has 
been bent so that its cross section forms an L-shape. 
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The 83rd iteration reconstructing two cracks. The collected data are the same as Figures 3a 
and 3b. The “true” crack locations are shown in dashed lines in the disk to the right. The 
currently active electrode locations are indicated on the circle to the left, 

6 Conclusion 

We believe that our experiments demonstrate the viability of Electrical Impedance Tomog- 
raphy as a method for locating defects in materials. The approach we have taken, using a 
priori information about the nature of the impedance profile, is clearly beneficial when such 
information is available. Moreover, our approach of using a few carefully chosen weighted av- 
erages of the boundary data, rather than all of the data as in least squares, further improves 
the performance of EIT. The adaptive movement of the electrodes makes the algorithm less 
likely to become trapped in locally optimal solutions which are not near the “true” solution, 
and since the electrodes converge to locations which make the Jacobian of the system far 
from singular, the algorithm should be less sensitive to noise in the voltage data. 
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The 6th iteration reconstructing one crack. The “true” crack locations are shown in dashed 
lines in the disk to the right. The currently active electrode locations are indicated on the 
circle to the left. 



The 31st iteration reconstructing two cracks. The collected data are the same as in Figure 4a. 
The “true” crack locations are shown in dashed lines in the disk to the right. The currently 
active electrode locations are indicated on the circle to the left. 
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The 5th iteration reconstructing one crack. The “true” crack locations are shown in dashed 
lines in the disk to the right. The active electrodes are indicated on the circle to the left. 



The 90th iteration reconstructing two cracks. The collected data are the same as in Figure 5a. 
Convergence has not been reached. The “true” crack locations are shown in dashed lines in 
the disk to the right. The currently active electrodes are indicated on the circle to the left. 
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Figure 6a 

The 6th iteration reconstructing one crack. The “true” crack locations are shown in dashed 
lines in the disk to the right. The currently active electrodes are indicated on the left disk. 




Figure 6b 


The 33rd iteration reconstructing two cracks. The collected data are the same as in Figure 6a. 
The “true” crack locations are shown in dashed lines in the disk to the right. The currently 
active electrode locations are indicated on the left disk. 



The 90th iteration reconstructing three cracks. The collected data are the same sis Figures 6a 
and 6b. Convergence has not been reached. The “true” crack locations are shown in dashed 
lines in the disk to the right. The currently active electrodes are indicated on the left circle. 



The 7th iteration reconstructing one crack. The “true” crack location is shown as a dashed 
curve in the disk to the right. The active electrodes are indicated on the circle to the left. 


24 











Acknowledgment s 


This work was performed while the third author was visiting the Universite Joseph Fourier 
during the academic year 92/93. He would like to thank his colleagues at “Laboratoire de 
Modelisation et Calcul” for the hospitality extended to him and in particular Jacques Blum 
for his efforts towards making this visit possible. 


References 

[1] Alessandrini, G., Stable determination of a crack from boundary measurements. Proc. 
Roy. Soc. Edinburgh, Section A 123 (1993), pp. 497-516. 

[2] Alessandrini, G., Personal communication. 

[3] Alessandrini, G., Beretta, E. and Vessella, S., Determining linear cracks by boundary 
measurements - Lipschitz stability. Preprint, 1993. 

[4] Andrieux, S. and Ben Abda, A., Identification de fissures planes par une donnee de bord 
unique: un procede direct de localisation et d’identification. C. R. Acad. Sci. Paris, Serie 
I, 315 (1992), pp. 1323-1328. 

[5] Barber, D. and Brown, B., Recent developments in applied potential tomography - APT, 
in Information Processing in Medical Imaging, Bacharach, S. ed., Nijhoff, Amsterdam, 
1986, pp. 106-121. 

[6] Bryan, K. and Vogelius, M., A uniqueness result concerning the identification of a 
collection of cracks from finitely many electrostatic boundary measurements. SIAM J. 
Math. Anal., 23 (1992), pp. 950-958. 

[7] Bryan, K. and Vogelius, M., A computational algorithm to determine crack locations 
from electrostatic boundary measurements. The case of multiple cracks. To appear in 
Int. J. Engng. Sci. 


26 



[8] Cheney, M., Isaacson, D., Newell, J., Simske, S. and Goble, J., NOSER: An algorithm 
for solving the inverse conductivity problem, Int. J. Imaging Systems and Tech., 22 
(1990), pp. 66-75. 

[9] Dia z Valenzuela, A. Unicita et stabilita per il problema inverso del crack perfettamente 
isolante. Thesis, University of Trieste, 1993. 

[10] Eggleston, M.R., Schwabe, R.J., Isaacson, D., Goble, J.C. and Coffin, L.F., Three- 
dimensional defect imaging with electric current computed tomography. GE Technical 
Report 91CRD039, Schenectady, NY. 

[1 1] Friedman, A. and Vogelius, M., Determining cracks by boundary measurements. Indiana 
Univ. Math. J., 38 (1989), pp. 527-556. 

[12] Gisser, D.G., Isaacson, D. and Newell, J.C., Electric current computed tomography and 
eigenvalues. SIAM J. Appl. Math., 50 (1990), pp. 1623-1634. 

[13] Liepa, V., Santosa, F. and Vogelius, M. Crack determination from boundary mea- 
surements - Reconstruction using experimental data. J. Nondestructive Evaluation, 12 
(1993), pp. 163-174. 

[14] More, J., The Levenberg-Marquardt algorithm: implementation and theory. Numerical 
Analysis (Edited by Watson, G.A.), pp. 105-116. Lecture Notes in Math. 630. Springer 
Verlag, 1977. 

[15] Nishimura, N., Regularized integral equations for crack shape determination problems. 
Inverse Problems in Engineering Sciences (Edited by Yamaguti, M. et al), pp. 59-65. 
Springer Verlag, 1991. 

[16] Nishimura, N. and Kobayashi, S., A boundary integral equation method for an inverse 
problem related to crack detection. Int. J. Num. Meth. Eng., 32 (1991), pp. 1371-1387. 

[17] Santosa, F. and Vogelius, M., A computational algorithm to determine cracks from 
electrostatic boundary measurements. Int. J. Eng. Sci. 29 (1991), pp. 917-937. 


27 



[18] Yorkey, T., Webster, J. and Tompkins, W., Comparing reconstruction algorithms for 
electrical impedance tomography. IEEE Trans. Biomedical Eng., BME-34 (1987), pp. 
843-852. 


28 




REPORT DOCUMENTATION PACE 


Form Approved 
OMB No 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 


1. AGENCY USE OUV( (Leave blank) I 2. REPORT DATE 

I December 1993 


3. REPORT TYPE AND DATES COVERED 


Contractor Report 


4. TITLE AND SUBTITLE 

RECONSTRUCTION OF MULTIPLE CRACKS FROM EXPERI- 
MENTAL ELECTROSTATIC BOUNDARY MEASUREMENTS 


6. AUTHOR(S) 

Kurt Bryan 
Valdis Liepa 
Michael Voeelius 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Institute for Computer Applications in Science 
and Engineering 

Mail Stop 132C, NASA Langley Research Center 
Hampton, VA 23681-0001 


S. FUNDING NUMBERS 

C NAS1-18605 
C NAS1-19480 
WU 505-90-52-01 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 

ICASE Report No, 93-96 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23681-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA CR-191584 
ICASE Report No. 93-96 


11. SUPPLEMENTARY NOTES 

Langley Technical Monitor: Michael F. Card 
Final Report 

Submitted to Proceedings of the 1993 SIAM Annual Meeting, Philadelphia, July 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

U nclassified-U nlimited 


12b. DISTRIBUTION COOE 


Subject Category 64 


13. ABSTRACT (Maximum 200 words) 

This paper describes an algorithm for recovering a collection of linear cracks in a homogeneous electrical conductor 
from boundary measurements of voltages induced by specified current fluxes. The technique is a variation of Newton’s 
method and is based on taking weighted averages of the boundary data. We also describe an apparatus that was 
constructed specifically for generating laboratory data on which to test the algorithm. We apply the algorithm to a 
number of different test cases and discuss the results. 


14. SUBJECT TERMS 

nondestructive testing; inverse problems; impendence tomography 


IS. NUMBER OF PAGES 

30 


16. PRICE CODE 

A03 


17. SECURITY CLASSIFICATION 18. SECURITY CLASSIFICATION 19. SECURITY CLASSIFICATION 20. LIMITATION 

OF REPORT OF THIS PAGE OF ABSTRACT OF ABSTRACT 

Unclassified Unclassified 


40-01-280- 


☆ U.S. GOVERNMENT PRINTING OFFICE: tfU - 528-064/86112 


tandar 

Prescribed by ANSI Std? Z39-18 
298-102 























